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The antiferromagnetic insulating cuprate Ba3Cu204Cl2 contains folded Cu02 chains 
with four magnetic copper ions {S — 1/2) per unit cell. An underlying multiorbital 
Hubbard model is formulated and the superexchange theory is developed to derive an 
effective spin Hamiltonian for this cuprate. The resulting spin Hamiltonian involves a 
Dzyaloshinskii-Moriya term and a more weak symmetric anisotropic exchange term be- 
sides the isotropic exchange interaction. The corresponding Dzyaloshinskii-Moriya vec- 
tors of each magnetic Cu-Cu bond in the chain reveal a well defined spatial order. Both, 
the superexchange theory and the complementary group theoretical consideration, lead 
to the same conclusion on the character of this order. The analysis of the ground-state 
magnetic properties of the derived model leads to the prediction of an additional non- 
collinear modulation of the antiferromagnetic structure. This weak antiferromagnetism 
is restricted to one of the Cu sublattices. 



PACS numbers: 75.50.Ee, 75.30.Et, 75.25.-l-z, 75.50.-y 



I. INTRODUCTION 

Electronic, superconducting and magnetic properties of the fast growing family of copper oxides and 
oxychlorides have attracted much attention last years. The undoped parent materials in the cuprate 
family are insulators showing a variety of low-dimensional magnetic properties. This variety extends 
from the quasi two-dimensional (2D) antiferromagnetic (AFM) behavioi§'i with possible admixture of 
a weak ferromagnetismE^ in the planar compounds to quasi-lD magnetic properties, observable in the 
chain-likeSi'i and spin- laddei0 systems. The superexchange theoryll3~lli provides the necessary basis 
which allows to derive and to estimate the main interactions, including anisotropic ones, responsible for 
the magnetic coupling of the copper spins. The form, and especially the magnitude of several interaction 
constants of the resulting spin Hamiltonian, however, strongly depend on pecularities of the Cu-O-Cu 
bond configuration in different cuprates (see Ref . |l^ and references therein) . One interesting class of mag- 
netic cuprates contains competing magnetic subsystems, like for instance tetragonal Ba2Cu304Cl2 built 
up of CU3O4 planes with two types of copper sites (Cul, CuII)0 It is well known that in Ba2Cu304Cl2 
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(as well as in Sr2Cu304Cl2) the moments of the Cul and Cull atoms order antiferromagnetically at 
different temperatures Tn.i w 330 . . . 380 K and Tnji « 30 ... 40 K, respectively.MB Below Tjv,/ a 
small spontaneous magnetization, Mq, within the basal plane of the tetragonal lattice has been reported. 
The corresponding magnetic structure has been analysed in Refs. ^-|7|, p^ , p^ . 

The aim of this paper is to develop the superexchange theory for the orthorhombic compound 
Ba3Cu204Cl2 (space group Pmma) with a Cu-O-Cu bond geometry rather unusual in the cuprate fam- 
ily, but having also two crystallographically different types of Cu sites.@ Actually, in this compound the 
edge-sharing CUO4 plaquettes are not aligned in the same plane, but form folded chains. The chain axis 
is parallel to the orthorhombic a-axis. The unit cell with two sorts, A and B, of crystallographically 
non-equivalent Cu sites is depicted in Fig. |^. Due to the typical Cu^+ (3d®) states the compound is an 
insulator. It was found to behave like a classical antiferromagnet with a Neel temperature Tjv of about 
Below Tjv and for applied magnetic fields Ha parallel to the a-axis this compound shows a spin-flop 
transition at /ioHa ~ 2.6T.0 Thus, the 'easy axis' of the antiferromagnetically ordered moments turns 
out to be the a-axis. Above the susceptibility follows a normal Curie-like behavior. From the Curie 
constant an effective paramagnetic moment of about 2 /is is derived, which is typical for Cu^+ in the 3(i® 
state. No weak ferromagnetism has been observedE Preliminary group theoretical analysis showedEl that 
in Ba3Cu204Cl2 weak ferromagnetism is forbidden while weak antiferromagnetism cannot be excluded. 
This symmetry analysis is extended below in close relation with the microscopical consideration based 
on the superexchange theory which yields Dzyaloshinskii-Moriya0@ type interactions. 

The organisation of the paper is the following. In the next Section the results of the bandstructure 
calculation in the local density approximation (LDA) for Ba3Cu204Cl2 will be briefly presented. The 
basic set of copper and oxygen orbitals responsible for superexchange processes will be fixed together 
with the estimates for the corresponding electronic transfer integrals and the orbital crystal field splitting 
parameters. At the next step, to provide a background for the superexchange theory the underlying 
multiorbital Hubbard model is formulated with spin-orbit coupling on copper ions involved. The per- 
turbation expansion of the multi-orbital Hubbard model is used in Section III to derive an effective spin 
Hamiltonian for the nearest-neighbour Cu-Cu magnetic interactions. A mean-field analysis of the model 
derived is also presented and a prediction on the ground-state magnetic structure in Ba2Cu204Cl2 is 
formulated. A complementary group theoretical analysis is presented in Section IV. Concluding remarks 
can be found in Section V. 
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II. BAND STRUCTURE AND THE UNDERLYING MULTI-ORBITAL HUBBARD MODEL 



The bandstructure of Ba3Cu204Cl2 was calculated within the local density approximation (LDA) using 
a recently developed full-potential nonorthogonal local-orbital minimum-basis scheme (FPLO)Ji Here, 
the Cu(3d, 4s, Ap), Ba(5s, bp, 5d, 6s, 6p), Cl(3s, 3p, 3d), and 0(2s, 2p, 3d) states were treated as valence 
states and the lower lying states as core states. For our purpose to extract tight-binding parameters 
it is sufRcient to perform a non-spinpolarized calculation. The non-magnetic solution shows metallic 
behavior with four bands crossing the Fermi surface, corresponding to the four copper atoms per unit 
cell. Hereafter the notations Bi, Ai, B2 and A2 will be used to specify the different Cu sites in a unit 
cell. 

Due to the folding of the Cu02 chains by nearly 90 degrees the bandwidth of the half-filled antibonding 
bands is relatively narrow (« 0.5 eV) and only half as large as in the planar (unfolded) edge-shared 
chains such as in Li2Cu02 or CuGeOs.B The orbital analysis for the four bands at Fermi level shows 
that these bands are built up mainly from the CnA-3dyz and the CuB-3da;2_y2 orbitals, respectively, 
with an admixture of all three 0-2px, -2py, -2pz orbitals (the global axes x, y, and z are parallel to the 
crystallographic a~, 6- and c-axes, respectively). The other Cu~3d orbitals give rise to a band structure 
(16 bands altogether) which spans the range of binding energies from 2 eV up to 4 eV below the Fermi 
level. The band states of predominantly O 2px, 2py, 2pz character fall into the range of binding energies 
from 2 eV up to 6 eV. 

The indirect coupling of the Cu-3c? orbitals via the intermediate 0-2p orbitals can be described by a 
set of effective hopping matrix elements. Considering first the top-most four bands formed by the active, 
i.e. spin carrying, CuA^3dyz and CvLB^3dr^2_y2 orbitals we obtained the following results by fitting those 
bands to a four-band tight-binding (TB) model.El The intra-chain hopping between the neighboring Cu^ 
and Gu-B sites, ~ 65 meV, is of the same order as the inter-chain diagonal Cusj^-Cusj hopping, 
t^^J^ ~ - 60 meV. Considering the hierarchy of transfer processes, the next terms were found to be 
t^^j^ ~ t^^g ~ 30 meV which correspond to the intcr-chain nearest neighbor hoppings along the y-direction. 
The inter-layer hopping parameters along z-direction are very weak. Therefore, the seemingly quasi-one 
dimensional compound Ba3Cu204Cl2 has to be classified, in a first approximation, as a two dimensional 
compound with spatially anisotropic interactions within a layer and weak couplings between the layers 
(which are necessary to explain the 3D magnetism and the experimentally found finite Neel temperature). 
That is true also for the exchange couplings of the spins related to the CvLA-Sd^z and GuB-3d^2_y2 basic 
orbitals. Below we will consider in great detail the microscopic origin of the dominant intra-chain nearest 
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neighbor interactions that involve rather strong magnetic anisotropics due to the non-trivial geometry of 
a particular chain. The remaining inter-chain couplings are expected to influence only the spin isotropic 
part of the exchange interaction since they occur in a more simple geometry. 

To develop the theory of superexchange, in the following the appropriate underlying Hubbard model 
is formulated for Ba3Cu204Cl2. Due to the S = 1/2 nature of the spins on the Cu+^ sites, a single-ion 
anisotropy cannot occur. A magnetic anisotropy can only be obtained by taking into account simultane- 
ously both, the spin-orbit coupling and the splitting of the orbitals by crystalline fields.0 The spin-orbit 
coupling of the copper ions is described by the -ffj^g term in the underlying electronic Hamiltonian. 
This term and the kinetic one, Ht, are considered as the perturbation while the zero-order Hamiltonian 
Hq — H^f^ + -ffo^'' is a sum of the on-site interactions at copper and oxygen atoms, respectively. Thus, 
the complete Hamiltonian is specified as follows 
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where the hole representation and the standard notation for cuprates are used. Here j, which is a 
composite index, denotes a cell number and the sort of the Cu site in the cell; ^ G j is used to denote 
an O site neighboring to the j-site. The term is the crystal field level for the m-th copper d-orbital 
at the j-th site. The results of the band structure calculationJH allow to estimate the crystal field 
splitting between the lowest e^Q (= £g) and the excited (m ^ 0) levels (~ ej^) as ef„ — Eq ~ 2eV 
= Ed- We neglect the difference of this splitting between the different excited m-states as well as between 
the nonequivalcnt Cu j-sites, which is of minor importance for the present purposes. For the further 
calculations only the copper on-site Coulomb integral, [/™™ , with m' = m is required, which is assumed 
in analogy with other cuprates to be [/™™ = Ud — S ^ lOeV. The crystal field splitting between the 
different oxygen |pifc)-orbitals (fc — x, y, z) located on a lattice site I is not taken into account in the 
present consideration. The estimate for the "bare" charge transfer gap Ap = — Eq ^ 4eV is deduced 
from the band structure calculations and a complementary finite cluster analysisH For the oxygen on-site 
Coulomb integrals Up'' = Up and Up'' = Up — 2Jp (for k ^ k'), the following estimates. Up ~ 4eV and 
Jp ~ 0.2 0.4eV, standard for the cuprates (see Ref. pSl and references therein), are assumed. 



Next, we consider the kinetie part, Hf . According to Slater and Koster (see Ref. ^) one can express the 
hopping amphtudes tjm,ik between the nearest neighbor copper \djm)^ and oxygen |p;fc)-orbitals as hnear 
combinations of the two parameters, (pda) and (pd?:). By using the approximation (pdn) ~ —l/2{pda) = 
—tp^/Vi which is vahd for transition metal oxides,0 one may write tjm.ik — Xmktp^- In this expression, 
for given \djm) and |p;fe), the factor Xmk is a function of direction cosines of the {j — I) vector. All these 
factors, of order unity, are calculated by using the routine procedure.S The only parameter remaining, 
tp^, is estimated from the bandstructure calculations. For the special composition of Cu04 plaquettes in 
Ba3Cu204Cl2 with only nearest neighbor Cu-0 hopping we obtained the following estimate t^^ ~ 0.5 eV 
by using the formula = (tp^)^/Ap. That value is smaller than the corresponding nearest neighbor 
Cu-0 hopping terms in other cuprateglll@0il since the effects of the different geometry, the neglect of 
direct 0-0 transfer and of the crystal field splitting at oxygen sites are all condensed into one effective 
parameter. The value A = 0.1 eV of the spin-orbit coupling, characteristic for other cuprates, is taken in 
the calculation. 

III. SPIN HAMILTONIAN AND THE MEAN-FIELD GROUND-STATE SPIN 

CONFIGURATION 

The derivation of the effective superexchange interactions using perturbation theory in the multi-band 
Hubbard model, has been presented already many times for planar cuprates.BEi'B Some peculiarities 
of the present derivation are due to the fact that the Cu02-chains in Ba3Cu204Cl2 are folded. The 
unit cell of this structure contains four different Cu ions which are denoted in Fig. 2 as Bi, Ai, B2 and 
A2. Therefore, there is a sequence of four different magnetic bonds (ij) = {Bi,Ai), (yli,i?2), (-82,^2) 
and {A2,B[), where the site B[ belongs to the neighboring unit cell. We start by deriving the spin 
Hamiltonian for a particular (AB)- or (BA)-bond of this sequence, and then extend the consideration 
to the entire structure. Two oxygen ions, {I = L,R), which are common neighbors of Cu^ and Cus, 
mediate the superexchange interaction between Gua and Cus. By applying an appropriate unitary 
transformation one can define molecular-type basic functions for the oxygen single-hole states in the form 
\pn) — {\pLk) ± \PRk))/ V2, with (n = 1, . . . 6 for k ^ x,y, z), that are used below. 

At the first step of the expansion procedure the excited |(ij>„)-orbitals (m ^ 0) are taken into account 
by additional vector hybrydization terms ^ Cj_„ • (t^^ connecting the ground state |djo)~orbitals with 
different |p„)-orbitals. Then, the original hopping processes ~ ijo.n and the additional new ones can be 
dealed simultaneously within the same approximation. The corresponding effective kinetic term is 



5 



H 



(AB) 



j=A,B n ap 



(2) 



where aa/3 represent the Pauh matrices and 
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Here, Lj,om are matrix elements of the orbital angular momentum operator and ij,„,„ are the original 
amplitudes for the hopping process between the excited m-th crystal-field d-state and the n-th oxygen 
p-state. 

Finally, the |p„)-orbital states can be eliminated in the perturbation procedure and the corresponding 
fourth-order processes are summed up into the following spin Hamiltonianll referring to one Cuyi-(O-O)- 
Cu-B magnetic bond 
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with the interaction constants 
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Here g„„' is given by 
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_Ud 2Ap 2Ap + U^^' 

where = Up and U™^ = Up — 2Jp (n' ^ n). Further on the notation (?„„ = g will be used. One 
may note a weak difference between gnn {n' — n) and gnn' {n' ^ n). Taking into account the particular 
form of the ijo.n {— tn,jo) and {— C* j) transfer parameters, we find the final expressions for the 
non-zero interaction constants: 



Jab = 4:9b AB - SJab, 



D\b = '89b ABC AB, 



(7) 
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where 
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1 - COS (f)AB) 



1 / A \ / pff\2 . , 
CAB = I — I [tpd J Sm (j)AB, 

6Jab « (■^) ^gble, (8) 



It is worth emphasizing at this stage that the vector Dab of the Dzyaloshinskii-Moriya interaction is 
directed along y-axis (i.e. — D\g — 0) and the symmetric anisotropy is described by only one non- 

zero parameter T^ab- Note, that the direction of Dab is in agreement with the rules given by Moriya.il 
The parameter (^ab is characteristic for the directed (yli?)-bond. More generally, a particular parameter 
(pij has to be assigned to each bond in the entire chain-like structure in the following way. Let us ascribe 
to the center of each j-th plaquette a local zj-axis perpendicular to the plane of the plaquette. For 
instance, in Fig. ^ the z{— zl) and z'{= Zj) are related to i — Bi and j = Ai, respectively. Then, 4>ij is 
the angle of the rotation which transforms the local z^-axis into the zj-axis. Therefore, one has 

'^BA = -^AB = -•/'Ia = ^ 84°, (9) 

where the left (right) indices refer to the left (right) atom in the directed bond. It can be directly checked 
that each magnetic bond is described by the same form of the spin Hamiltonian (^-(||). By using the 
relations (0-|l) it can be easily seen that the parameters J and F^^ are bond-independent while the 
component D^j changes the sign in accordance with (^. The spin Hamiltonian for the entire chain-like 
magnetic system can be written in the following form 

= 4g {9S^ ■ Sj + 2bc£,^jd- S^ x Sj + {2(d5,)(d5j) - S, ■ 5^}) , 

sjiiv) ^ -sjsjj - sryy \^2{dSi){ds^) - sj^} , (lo) 

where d is the unit vector along the global fe-axis. ^tj = sgn^^j determines the space pattern of the 
Dzyaloshinskii-Moriya vectors for the chain-like structure (see Eq. (9)). 

Below the arguments given by Shekhtman, Entin-Wohlman and Aharonj^Il for the single-bond su- 
perexchange interactions are extended to the chain-like magnetic system with this special pattern of the 
Dij vectors. Actually, the spins in the lattice can be subdivided into two subsystems in such a way that 



the first (second) subsystem is formed by the spins in B (A) positions only. Let us now introduce the 
following redefinition of the spin variables 
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— k sin 9 
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X d 
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+ k sin 9 


Sa2 


X d 



(11) 



where k = ±1 and ktaii9 = — 26c/(6 — c ). It should be noted, that the transformation ( |ll| ) corresponds 
to a rotation of the spins Sai iSA2) by an angle 9 (—9) with d as the rotation axis. The Hamiltonian 
Hq''^ is strictly transformed into an isotropic formil while concerning the remaining STi^^^^ an additional 
justification should be made. Actually, no terms higher than of second order in (A/e^) have been kept 
up to now. By noting that 5J,5T ^ (X/edf' and in ( |ll| ) the angle 9 ~ (A/sd), the consistent way is to 
keep the form of STi^-^^^ but using the transformed spin variables. Finally the transformed Hamiltonian 
n = J2 ^^'^^ takes the form 



H^*^) = JS^SJ ~ f \^2{dS^){dSj) - (^S^SJ^ } 



(12) 



with the renormalized constants J = 49(6^ + c^) — 5J, and f = ST^y > 0. 

A mean field analysis of this Hamiltonian shows that the expected classical ground-state configuration 
is a collinear antiferromagnetic array of the transformed spins (|l^) in the magnetically easy xz-plane, 
f > 0. We emphasize that according to ( [l^ ) the staggered moment is not confined to a particular 
direction in the easy a;z-plane. This 'residual' symmetry of the derived superexchange can be broken by 
additional interaction terms not included into the consideration up to now. Postponing a discussion of 
possible sources for additional anisotropics let us now assume that the 'residual' symmetry is broken in 
such a way that the staggered magnetization is parallel to the x-axis (i.e. crystallographic a-direction). 
It means that in the ground state the spins Sbi {— Sbi) and Sb2 {— 'SB2) are aligned in the af-direction. 
In the mean-field approximation one may immediately write 



Sbi — ^Sai — Sb2 — ^Sa2 — Sex 
and performing the inverse transformation one obtains for the original spins 

Sbi — Sb2 — Stxj 

Sai = ~S [cos 9 ■ ex + k sin 9 ■ Bz], 

Sa2 — ~S [cos 9 ■ ex — k sin • e^] . 



(13) 



(14) 
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The corresponding picture for the double degenerate classical ground-state spin configuration is given 
in Fig. 3. A weak transverse (along z-axis) modulation in the second sublattice is imposed on the 
strong antifcrromagnetic correlations between the spins belonging to different sublattices. This weak 
antiferromagnetic modulation is entirely due to the Dzyaloshinskii-Moriya interaction. The angle 9 for 
this modulation is estimated to be 

with an absolute value of \Dij \ ^ 0.25 meV corresponding to J ~ 5 meV. It is worth mentioning that the 
magnetization measurements! for Ba3Cu204Cl2 have clearly shown that the crystallographic a-direction 
(x-axis in our notations) is the preferred direction for the staggered moment below Tn, where Tjv — 20K is 
the temperature of the three dimensional antiferromagnetic ordering in Ba3Cu204Cl2. According to this 
analysis it is expected that the ground state in this compound is not a simple coUinear antiferromagnetic 
configuration, but involves also a weak antiferromagnetic superstructure due to the Dzyaloshinskii-Moriya 
interaction. 

To our knowledge, in cuprates only weak ferromagnetism (WFM) was reported up to now. The 
corresponding spin canting angle, which measures the deviation from coUinearity, is rather small, ^'VVFM — 
0.05° (compare for instance Refs. p|, p6p^ , |3^ ). In the present case the spin canting angle (p^), i.e. ~ 3°, 
is more than an order of magnitude larger. The reason for this difference is the following. Although the 
spin-orbit coupling (\/ed) ^ 0.1 is nearly the same in both cases, the geometrical factor is much larger 
in the present case as |sin(/),y| ~ 1 due to the strong folding of the chains {\(t>ij \ — 84°). 

It is necessary to discuss the possible infiuence of transverse inter-chain magnetic interactions to show 
the validity of the derived spin anisotropy. We assume that the dominant interaction between two 
neighboring spins belonging to different chains in a layer is the isotropic superexchange with an AFM 
interaction constant J^. By using the results of the bandstructure calculations (Section II) we found the 
estimate J_l/ J ^ 0.25, where J is the intra-chain exchange constant. Due to the presence of the center 
of inversion for the transverse Cu-Cu magnetic bond a Dzyaloshinskii-Moriya term cannot occur and 
the symmetric anisotropy is expected to be very weak. The bandstructure results indicate also a weak 
inter-layer exchange Jz ^ J_\_, J and more distant inter-chain Cu-Cu exchange interactions (for instance 
to third neighbors in the plane of Cu^). However, the spin anisotropy of the latter interaction is expected 
to be small due to the planar geometry of the corresponding exchange path. 
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IV. SYMMETRY AND SPATIAL ORDERING OF THE DZYALOSHINSKII-MORIYA 

VECTORS 

In order to derive the spatial order of the Dzyaloshinskii-Moriya vectors by symmetry considerations, 
at first the symmetry analysis of the antiferromagnetic states in Ba3Cu2 04Cl2 (following Refs. p6[-p8|), 
presented elsewher^^, will be summarized. It is assumed that the lattice constant a of the magnetic 
unit cell and that of the crystallographic unit cell are identical. As has been pointed out in Rcf. ^ it 
is sufficient to consider the independent symmetry operations. For the space group Pmma besides the 
identity there are only three independent symmetry elements. 

A very convenient choice for these elements is: rotation axis C2y through the point (xo=0,zo=l/2), 
rotation axis C2z through the point (xo=l/4,yo=0) and inversion I located at (xo=0,yo=0,zo=l/2). As 
these elements map a given chain onto itself no assumptions concerning the relative orientation of moments 
in different chains are necessary. The magnetic moments reside on the 2c-sites Cvlai, Gua2 and Cuy^'i, 
and on the 2e-sites Cusi, Cvlb2 and Cus'i (see Fig. ^ in which also the mentioned symmetry elements 
are sketched). Obviously the ion Cu^n (Cusi) is equivalent to Cu^'i (Cus'i). Because the sites of the 
2c- and 2e-sublattices are crystallographically not equivalent, the magnetic order has to be characterized 
by two antiferromagnetic (AFM) vectors 

La = Mai - Ma2 Lb = Mbi - Mb2 (16) 
and two sublattice magnetizations 

ruA = {Mai + Ma2)/2 tub - [Mbi + Mb2)/2 (17) 

where Mai {Mbj) are the magnetic moments of the Cu Ai {CuBj) ions at the sites A^ (Bj). The interchange 
of the moments Mi <-> Mj generated by the symmetry operations is given in the second and third 
column of Table I. As only one ordering temperature is found, both AFM vectors have to transform 
according to the same magnetic group, compatible with the crystallographic space group Pmma. The 
transformation behavior for each component La.x, LA.y etc. and for the sublattice magnetizations rfiA and 
ms is described in Table I by two numbers characterizing two steps of the corresponding transformation: 
the transformation of La, Lb etc. due to Mi Mj (first number) and the remaining part of the 
transformation (rotation, inversion; second). The entry +1 means no change of the considered vector 
component, -1 means a change of sign. In the line below the resulting magnetic symmetry element is 
given. As I La = La, but RIL^ = Lb the AFM vectors L2c and L2e cannot simultaneously become 
different from zero. As each component of rriA and rfiB transforms according to the same magnetic group. 
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the only possibility for antiferromagnetic ordering is — toa.k - tt^b,x Ly ^ '^A,y - ^B,y ^ Q ov Lz 
= fTT-A.z - TTT-B.z 7^ 0. This means, ferromagnetic order is predicted for each of the two crystallographic Cu 
sublattices. The moments of the crystallographically different sublattices are antiparallel to each other. As 
different components of rfiA and rfiB belong to different magnetic groups (Table I) (i) only one component 
of the AFM vector L may be different from zero, i.e. L is parallel to one of the crystallographic axes, and 
(ii) weak ferromagnetism is excluded, contrary to the case of Ba2Cu304Cl2. This is in agreement with the 
experiments. The symmetry analysis cannot predict the direction of the AFM vector. Experiments show 
Lx =/= 0.0 For that case, independent magnetic symmetry elements are E, RC2z, RC2J, and 1. Interestingly, 
Table I shows, that for that case a small La,z 7^ i.e. weak antiferromagnetisrrS can additionally occur. 

As shown in Sect. Ill, within the unit cell this weak antiferromagnetism is described by a sum of 
Dzyaloshinskii-Moriya terms 

HwAFM — D^j^ ■ {Mbi,xMai^z — Mbi^zMai^x) 

+ -^AB ■ {Ma1,xMb2,z — Mai,zMb2,x) 
+ -DfiA • {Mb2,xMa2,z — Mb2,zMa2,x) 

+ D^2b' ■ {Ma2,xMb'1,z - Ma2.zMb'1,x)- (18) 

This interaction has to be invariant with respect to the symmetry operations RC22, RC2J, and I. Using 
the transformation properties of the moments according to Table I, the equivalence of sites differing by 
a lattice translation along the a-axis, results in 

Dh\ = -D\% = -Dl\ D^2b' ■ (19) 

As to be expected, the spatial order of the local Dzyaloshinskii-Moriya vectors determined from the 
symmetry of the magnetic unit cell is the same as that derived by analyzing the relations between the 
different bonds in the previous Section. Moreover the symmetry analysis together with the experimental 
results shows, that in the transformation described in Sect. Ill (canting of the moments) the GwBi 
moments have to be fixed, as Lb has to remain zero for an AFM state with moments parallel to the 
a-axis. 

V. CONCLUSIONS 

The superexchange interaction between two neighboring but crystallographically non-equivalent Cu 
ions in Ba3Cu2 04Cl2 was analyzed within a multi-orbital Hubbard model. This superexchange inter- 
action was expressed by an effective spin-spin Hamiltonian. As there is no center of inversion for the 
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considered Cu-(0-0)-Cu entity, this interaction involves the Dzyaloshinskii-Moriya term, which refers to 
the corresponding bond. By means of an analysis of the geometric relations between neighboring bonds 
(plaquettes) the spatial order of these locally defined Dzyaloshinskii-Moriya vectors was determined. This 
spatial order was also derived by a symmetry analysis taking into account the experimentally found 'easy 
axis' of the staggered magnetization. These considerations revealed, that spin canting occurs in the sub- 
system of the Cua magnetic moments only. Within the presented microscopic model, for the magnetic 
moments all directions within the a-c-plane are equivalent. The anisotropy, leading to the experimentally 
observed spin-flop transition for applied fields parallel to the a-axis could not yet be explained. Further 
investigations have to be done to find the microscopic origin for additional anisotropics, which break 
the easy-plane symmetry of the superexchange model and explain the experimentally observed easy-axis 
behavior. With this respect, two main contributions have been ignored in the considered Hamiltonian 
(|l]). First of all, the theory of superexchange can be developed at a more sophisticated level by adopting 
more details of the electronic structure of the real compound. For instance, the form of the symmetric 
anisotropy tensor is rather sensitive to an actual crystal-field splitting of p-orbitals on oxygen ions. This 
splitting should be taken into account into the theoretical scheme. A second reason leading to a break- 
down of the easy-plane symmetry involves the direct-exchange contributionil to the symmetric anisotropy 
term. This is due to the exchange part of the two-site Coulomb multi-orbital d ~ d correlations which 
should be also incorporated into the Hamiltonian (^. However, for such calculations reliable quantita- 
tive estimates for the parameters are required (i.e. the crystal-field states of the 0-ions and the two-site 
exchange Coulomb integrals of Cu-ions). Since the direct exchange contribution can be hardly obtained 
within the LDA band structure analysis, we have to leave this problem for future investigations. From 
an experimentally point of view, the detection of the discussed weak antiferromagnetism is a challenging 
task. 
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FIG. 1. Two unit cells of Ba3Cu204Cl2, space group Pmma. The CU-O4 plaquettes form folded chains with 
their axes parallel to the o-axis. To better demonstrate one of these Cu02-chains, the origin in the figure was 
taken at (0,1/2,0). The crystal structure as well as the chains contain two types of copper sites, Cua and Cus. 

FIG. 2. Spatial ordering of the Dzyaloshinskii-Moriya vectors -D^^ etc. related to a pair of Cu atoms located 
at Bi and Aj. (The first indices refer to the left atom in the pair). The z~ {z'-) axes are locally defined and 
are perpendicular to the corresponding plaquette. The dashed line is the local z- (z'-) axis of the plaquette left 
to that under consideration and is used for the definition of the angles 0^^: The rotation of the corresponding 
dashed line onto the local axis results in the rotation angle i^^^ , the sign of which defines the direction of the local 
Dzyaloshinskii-Moriya vectors 5^^. The crystallographic axes a, b, and c form the global coordinate system. 

FIG. ,3. Within a classical picture (mean field approximation) the Dzyaloshinskii-Moriya interaction between 
the magnetic moments (arrows) of the Cua and Cus atoms results in a small canting of the Cua moments (canting 
angle ±9). The ground state is two-fold degenerated (a,b) and reveals weak antiferromagnetism. The two states 
differ by the sign of the canting angle and the direction of the of the additional antiferromagnetic vector La,z = 
Mai,z - Ma2,z of the Cua sublattice. x and z are global axes equivalent to the crystallographic a- and c-axes. 

FIG. 4. Symmetry elements of a CuA-Cus-chain in Ba3Cu204Cl2. The magnetic moments are located at the 

2c-sites (Ai, A2 and A'l) and at the 2e-sites (Bi, B2 and B'l). The symmetry elements mapping the chain onto 
itself are: axes C2z going through the B-sites, axes C2J, (sketched by small ellipses) perpendicular to the a;-«-plane 
and going through the A-sites, and inversion centers located at the sites Aj. 
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TABLE I. The first column lists the independent symmetry elements of the space group Pmma, here chosen 
to be E, C2z, C2y and I. The second and third columns give the magnetic moments Max and Mbj (i, j = 1, 2) 
generated from Mai and Mbi by applying these symmetry elements. In the remaining columns for each component 
of the antifcrromagnetic vectors La ~ Mai — Ma2, Lb = Mbi — Mb2, and the two sublattice magnetizations in a 
= {Mai + A/a2)/2 and ms — {Mbi +Mb2)/2 two numbers are given: the first takes into account the interchange 
Mai «-> Ma2, Mbi <-+ Mb2 (1 <-» no change of the vector components, -1 <-> vector components change sign). The 
second +1 or -1 describes the remaining part of the considered transfomation (rotation, inversion). The resulting 
magnetic group elements are given in the line below these two numbers (R - time reversal). 
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